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I.  HmWDUCTIO* 


¥hit  report  lUKuiriMi  tha  praaant  atatua  of  a  eontinulnit  Inraatlgatloo 
Into  tha  fQfstariaa  of  tha  sMll-seala  atnoapharle  rortax*  uaually  knovn  aa 
tba  "dust  da^I"  or  "vhirlvind."  Nny  of  tha  baalo  coneapta  of  Tortax  fJ.ov 
praaontad  balov  appX^  to  tha  largar-aeala  toropdo  aa  vail.  Tha  rafaranca 
litaratura  on  tha  theory  of  anall-eeala  rortleaa  In  aataorologl eal  Journals 
la  relntivalj'  ai'ama}  bouavari  thara  ia  an  axtanslTr  litaratura  on  basic 
vortex  BCtiona  available  In  Journals  on  fluid  aMchanles. 

Tho  first  application  cf  tha  ravlar-^tokas  aquations  of  fluid  sotion 
to  a  viscous  vortex  problen  saavs  to  have  bean  tha  invastigationa  of 
Taylor  (1918)  and  Terasawa  (1922),  suiraarizad  in  Drydan  (19S6)i  on  tha 
decay  of  a  vortax  consisting  only  of  tangential  valooitiaa.  Bursars  (19^^ 
19^6),  Rott  (1958-1959)*  Sullivan  (1959)*  Donaldson  and  Sullivan  (1980)* 

Long  (1958-1961),  Lavallan  (1962)  and  Wbbb  (1962)  have  invaatlgatad  etaady- 
stata  vortax  flovs  that  are  axisynratric.  Cf  tbasa,  Lavallan* a  traataant 
is  tha  noat  coorrahensiva.  Outaan  (195T)  and  Kuo  (I96H}  have  worked  out 
aolutioos  vhlch  incorporate  the  affects  of  tbamal  energy  input  on  tha  flow. 
Boundary  layer  affects  have  bean  dlaeusaad  by  Rott  (1962),  Mataorolof.ical 
appllcetion  has  bean  nada  by  VauAhn  (1928),  Riaphreya  (19^0),  Villiana  (19^8), 
Battan  (1958),  as  vail  as  KUo  and  Vabb.  Sinclair  (I96h)  has  uada  sosm 
Masuraaents  of  tanparatura  end  velocity  profiles  in  dust  devils  using  a 
portable  instruBant  probe. 

The  instability  of  rotating  flovs  has  bean  invastigatad  by  a  ntadMr 
of  people,  including  Harrlaoo  (1921),  Taj'lor  (1983)*  Synra  (1938),  Lin  (]9>5) 
end  Chandrasddiar  (1962);  hut  their  at'^as  have  dealt  prlnarily  vlth 
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Couitt*  flov.  So  for  M  th«  otithora  kxKW,  no  eotq>r«honslr«  iQr»stlcatlc& 
of  vortox  stability  Itidspsndant  of  bounding;  vails  has  been  carrlsd  out. 

la  this  report,  aone  significant  raiults  of  thtsa  prior  Inrestlga- 
tlona  era  briefly  presented  and  dlacuased,  and  vhere  appropriate  rrer»erallsa- 
tioDS  of  these  results  are  offered.  A  partial  analysis  of  the  stability 
of  vortex  flov  is  also  included.  Mnall;',  e  prograai  of  future  Invest igac 
tloas  Is  outlined. 


II.  STEADY-STATE  VERTEX  FLOU 

The  starting  point  for  all  dlscuLilcna  of  lanlnar  viscous  vortsx 
flw  In  the  set  of  squatloas  of  fluid  r.otlon  (Uavisr-Ftoksa)  expressed 
in  cylicdrloal  ceordlnatss,  e.g.  Hinse  (IPS 9)  PP  22-23: 
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Here  r,  S  and  z  are  the  radial,  angular  and  axial  coordinates 
respectively;  t  is  tbe  tliae;  u,  v  and  v  are  the  conponents  of  velocity  in 
the  r,  0  and  i  directions  respectively;  and  p  is  the  pressure, e  the  Oen.lty, 
V  the  kineaatlc  noleeular  viscosity  and  g  t.ie  acceleration  due  to  cr 
The  operators  D  and  ^  ars  dsfined  as 


Di'u|-.S»  .v|- 
Sr  r  *6  9t 
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and 


2  ,*.2 
r  36 


In  the  foiJoving  diPcuriioDa »  etce^t  ^or  the  study  of  t^ell  fer-tui'bi^ 
tinne,  the  nct^cn  it  aseuand  to  be  exlsi-metrlc,  so  that  the  derlvstlvei 
vith  resijec^  to  0  vealsh  In  the  nbcve  equations,  resulting  in  considerable 
slmFllfl cation,  Purtherr.oi^,  In  a  first  approach  to  understandinr  the 
Bto^.ions  In  a  ** dust-devil,"  it  is  convenient  to  eosuae  that  the  horiiontal 
t'eioclty  contponents ,  u  end  v,  dc  not  varj'  Ir-  the  vertical.  Thla  ne.^lecta 
the  problem  of  flcr.'  in  t!:e  curface  boundar:.''  layer,  but  ahould  be  appro¬ 
priate  for  the  flow  pattern  at  some  distance  abc/e  the  surface.  Eouationa 
(2.1),  (2.2)  and  (2.3),  usinp  th?  above  eeevnptlor.s,  besone 


where  u  ^  u(r),  v  ■  v\r),  w  »  w(r,s),  p  »  p^r,a),  v  ■  conatant  end  o  ■o(») 
constant.  It  ie  ass'iF^d  that  the  earth 'a  rotation  has  little  effect  on 
the  size  of  disturbance  beinr  considered,  hence  the  Coriolis  terras  hv^o  boon 
neglected.  It  is  further  aaruned  that  the  viscosity  is  constant:  .his 
eecvrption  will  be  adsquste  fcr  rolecular  viscoeity,  but  is  in  dou'.  ^  vhea 
eddy  viscoeity  is  considered. 

The  eQuatiOQ  of  nase  contlnxiity,  assuming  that  the  density  is  approri- 
oatelv  constant,  may  be  written  ae 

^  1^  (ru)  ■  -  »c(r)  2  -  IJ,  or  u  ■  -  p  /  r*  (r)  dr 

0 


(2.T) 
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vt)«rtK(^)  is  ths  horlsontsl  conrercrsncs ,  s  functloD  of  ra^us  only.  With 
ths  IntroduetioD  of  this  convsrffsncs  tsra,  1>*«  flow  omy  be  coosldered  tvo- 
dlssnsional  in  the  r,§  piano >  •t(r)  then  plnye  the  role  of  a  distributed 
velocity  sish* 

The  vertical  cooipooent  of  vertlflty,  the  only  cocpor.ent  present  in 
this  tvo<-dlDsn8ioQal  setup,  is.  In  cyllnlrieal  coordinates, 

C  ■  ~  (rv)  (2.8) 

Substitution  of  (2.6)  in  (2.5)  rives 

(2.9) 

as  noted  by  Laab  (1932,  p.  579)*  Integration  yields 

C  "  Co  exp  o/^  HP"  *  central  vorticlty  (2.10) 

and  a  second  Integration  yields 

rv  ■  4o  cj^  **  drj  dr  (2.11) 

For  various  special  choicos  of  u(r),  equation  (2.11)  oey  be  solved 
aoalytlcally  for  rv  and  v.  One  such  case  is  when 

u  ©  (a.r:) 

Then  (2.7)  gives  k  (r)  ■  <©  ■  constant,  stating  that  horliontal  con'"^«riPCiice 
is  ^J&lform  for  this  case.  Such  an  assumption  nay  hold  quite  well  over  a 
large  region  near  the  central  axis  of  the  vorter.  The  resulting  •’■•iiT'icDtifll 
velocity  profile  Is 


where  is  the  circulation  at  r  •>  *•.  This  is  essentially  the 

solution  of  ftirgers  (19^0)  end  Bott  (1956).  If  v  is  defined  as  the 

ID 
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maxinua  tansantiiLl  wloeity  sad  r  th«  corraipocdlng  radius  at  vhleh  tha 

Q 

naxlmun  ocoura,  tban  dlffartDtiation  of  (2.13)  yialda 

2 


•  •^-.2.51 


(2.\K) 


as  a  charactarlatlc  autibar  for  a  vortex  ^'Ith  this  particular  velocity 
profile t  as  implicitly  stated  by  Rott.  rXirwharnore, 


.72. 


r  5  2trr  r 

B  IQ  B 


(2.15) 


a  ♦  1 

indicatias  tha  lavlatloc  of  this  velocity  profile  fron  the  simple  con¬ 
servation  of  angular  moiuentun,  F  ■  constant,  one.  If  nolecolar  viscosity 
of  air,  O.lU  cn^  sec  is  used  in  (2.1^),  then,  for  a  typical  dust- 
devil  dimension  of  r  ■  100  cr,  the  corvergenca  kq  Is  found  to  be 

In 

7  X  10“^  sac  If  the  convargin?  air  starts  at  a  radius  of  r.  and 
moves  Invard  to  a  radius  in  time  t,  this  radial  inflow,  obtained  by 


integrating  (2,12)  fron  r^^  to  r,,  is  given  b!' 


r-/r-  ■  or  t  ■~lcl(r,/r„)  (2.l6) 

c  l  *  rjg  i  *: 

I* 

For  ^2.^*2  "  ^  ^  ^  about  5  hours,  vhlch  is  unrealistic. 

If  tha  acleculr.r  viscosity  is  raolaced  by  an  effective  eddy  vlseoalty  of, 
e.g. ,  10^  cn^  sec"^,  "  5  x  lO"^  aec**^  for  r^^  "  100  cm,  and  for 
t^/rg  ■  2,  t  »  28  sec,  c  ror.sonsole  period. 

The  eddy  viscosity  should  be  taken  et  least  in  the  outer  portion  of 
the  dust  devil  but  probably  not  in  the  center  part  wbtT-e  solid  rotation  is 
approximated.  Pott  argues  alRllarly  for  the  application  to  a  toruado. 

The  vorticity  for  the  above  cese  asaunee  a  Oaussian  distribution: 


C  •  Co 


-iCor^/l*u 
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* 


6- 


Anoth«r  eaM  of  latoroot  li  that  for  «ihleh  tb*  coovtereenc*  1b  zero  out¬ 
side  of  a  "sink"  ref  loo  of  radius  r^^.  For  this  case  u  "-b/r  for  r  >  r^^,  and 
the  resulting  angular  monentun  distribution  is 


rr 


(2.18) 


vhlch  Is  tba  Soluti<Hi  obtained  by  Webb  (19^)  for  the  ssiae  sink  assumption. 
Vor  rv  to  be  a  bounded  function  for  large  r,  it  is  reouired  that  b/v  *2, 

Fbr  b/v  *2,  the  invlscld,  solution  rr  ■  results.  The  ratio  of  the  cir¬ 


culations  r_  and  r. 


r  -  -  2 

:i  _  V  ^ 


2eriVi  is 


r*b 


V 


-  1 


The  beharlor  of  the  flow  vlthln  the  sink  region  must  be  specified  by 
eosM  other  condition,  such  as  lulform  convergence,  etc* 

A  third  analytic  case  of  some  interest  is  one  in  which  u  is  a  constant, 
nils  unrealistically  implies  infinite  convergence  ai  d  vertical  acceleration 
alcog  the  central  axis,  but  may  serve  as  a  limiting  case  for  vortices  where 
the  convergence  and  updraft  are  highly  concentrated  near  the  center,  in 
contrast  to  the  uniform  convergence  case.  With  u  ■  constant,  the  resulting 
velocity  profile  is 

T  -  ^  [i  -  (x  *  ssj  r.  .  (?.ao) 

and  where,  in  the  same  manner  as  in  the  first  solutlcn  and  (2.13), 


ur_ 


8.  J,  1.80, 


*  B  *  1 


0.5l» 


(2.21) 


The  velocity  profiles  d^ven  by  (2.13)  and  (2.20)  have  very  slnller  shapes 
despite  the  great  contrast  in  the  respective  convergence  patterns,  see  figure 


Of  cour8«i  th«i*  Telocity  prefix  ••  say  be  pieced  together  (  for  exeaple. 


u  «  for  O'  ^  r  <  P,  u 


for  r.>  r,  then 


rr  • 


(2.22) 


for  r  >  P,  vhereas  equatloQ  (2.13)  applies  for  r  <  R.  Ultiaiately,  the  only 
Justification  for  constructing  such  profiles  is  that  they  are  readily 
obtainable  in  analytic  fom,  vhich  does  not  ing>ly  that  they  are  veiy 
good  approxlBatloas  to  reality. 

A  flow  pattern  which  is  sosMr^at  nore  general  than  the  above »  and 
which  presusably  has  greater  correspondence  with  reality,  is  one  In  which 
the  horiaontal  canvergence  is  assumed  to  be  a  ’^lll-sbaped"  function, 

^(r)  •  -r—Bg-  (2.23) 

a*^  +  r 

where  a  is  the  radius  at  vhich  k  falls  to  Ko/2.  The  correspoodinp  function 
for  radial  velocity  is  then 


(2.2U) 


The  maximum  radial  velocity  is  readily  shown  to  occur  at  r  ■  1.996  a  and  has 

the  value  u  *  0.U02  Koa.  The  tangential  velocity  cannot  be  obtained  in 
m 

analytical  form  for  this  ease,  but  it  can  readily  be  computed  numerically 
using  finite-difference  Integration  methods.  For  purposes  of  conpc;’:* non , 
the  following  dimensionless  quantities  are  defined.: 


•8* 


KoTji 

X  t  r/r^,  j  5  t/v^,  fi  1  a  i  — ^  (2.25) 

vtier*  Y.  Ifl  tbt  TMXiaMD  tangautial  Ytloelty  and  r  tb«  corr«*pocdlng  radius. 
Iha  paraiaetar  <  glYes  a  aeasura  of  tba  eoceantretlon  of  eoDYorganca  naar  tha 
eantral  axla;  for  vary  aaiall  6  tha  eonvarganea  la  highly  eoncaotratad  along 
tha  aadat  for  <  *•'  «•  tha  eonvargaBoa  profile  approaehaa  unlfomity.  Baaulta 
of  tha  nitaarleal  caleulatlona  for  c/iCg  and  y  aa  functiona  of  x  are  plotted 
In  figure  (2).  Ihe  relation  hetvaen  a  and  6  la  plotted  in  figure  (3), 
ahovlng  that  a  2.51  aa  ^  It  la  noteworthy  that  for  values  of 

^  *  It  the  deviation  of  the  tangential  velocity  fron  the  uniform  convergence 
aolutloB  (2.13)  la  vary  small. 

Nora  generally,  tha  eonvergenea  can  ha  expressed  in  the  fora  of  a  aeriea 

n  a,  icq, 

•'(r)  «  I  g— -o- 

1-1  af  ♦  r*' 

vbara  tha  a^  and  are  individually  different  values  of  half-vldth  and 
central  convergence.  For  n  ■  2,  a^^  <  Sg,  k©!  ^  0  and  Xog  >  0,  a  concentric 

two-^all  vortex  results,  vlth  divergence  dominant  in  a  central  "aye"  region 
and  eonvergenea  elaevbere.  Such  a  profile  Is  possibly  characteristic  of  the 
tornado.  Gutman  (1957),  Sullivan  (1959),  Donaldson  and  Sullivan  (i960)  and 
Kao  (1964}  have  obtained  solutions  for  concentric  tvo-cell  vortices.  For  all 
tha  foregoing  solutions,  the  pressure  profile  is  c^lven  very  closely  by 

r  2 

p  -  Po  ♦  p|  (^)df  (2.27) 

ebera  Is  tha  central  pressure,  p  tha  density  for  the  heicht  of  Interest 
and  v(r)  the  tangential  velocity  wofila  for  tha  particular  vortex.  The  othsr 
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t9rm  of  oauAtloB  (2<lt)  upoB  Igto^otloa  ftw  •  lEt^l  corr«ctioB  that  rofulta 
in  roduelDg  the  pressure  gradient  given  by  (2.27)  slichtly. 

The  diasipation  of  energy  in  the  vortex  due  to  viscous  forces  has  been 
discussed  by  Burgers  (19^8)  and  Rott  (19?9)*  Tott  gives  the  energy  eqviation 
for  axisymaetric  flov  in  the  foUovlTig  foiri: 


PUT  g  .  »v»  ♦  i  fj.  (tr  g)  (2.28) 

where  S  is  entrophy,  T  terperature,  k  thermal  conductivity,  and  •  the  viscous 
dissipation  function  given  by  (Lanb,  o£.  clt»»  p.  580) 


*  •  2(^)'  -  >  2(U)'  *  (||)’  ♦  (g  *  |S)%  (r  ^  (I)r(2.29) 


•Jw 


For  Burgers*  uniform  convergence  solution,  the  resulting  values  are 


«  « 


rSl  - 1) 


♦  6k. 


(2.31) 


px  r  2 

D  ■  -j—  (2.32) 

as  given  by  Burgers.  The  approximate  equality  sign  indicates  here  that  the 

2 

6kq  term  is  negligible  in  the  integration  of  *  cut  to  a  radius  much  larger  than 

r  .  Thus  the  dlsEipation  function  is  primarily  a  function  of  the  tangential 
m 

shear  (dv/dr  -  v/r)  and  the  res'jd.ting  dissipation  depends  on  the  **basie 
circulation,"  as  well  as  the  coa%'er6ence.  It  can  be  shown  that  the  dominant 
dissipation  term  is  always  the  tangential  shear,  and  that  (2.32)  is  a  good 
approximation  for  the  other  vortices  discussed  above  as  well. 

Equation  (2.28)  can  now  be  used  to  calculate  the  (Synamlcally  irdured 
temperature  difference  between  the  center  of  the  vortex  and  the  enviru  - 
ment,  as  is  done  by  Eott  for  a  perfect  gas  with  various  thermal  condu.-tlvit'.' e . 
>'or  the  limitiBg  ease  of  no  conductivity,  be  finds 


b 


10- 


T_  -  To  ■  ^  -R  -g  (2.33) 

®p  I6ir^r  ^ 

«tMr«  To  it  tilt  ctntml  totptrature  and  la  the  ttoparatura  at  r 
Htnet  tht  motion  rtaulta  in  eoolins  at  tha  air  approachea  tha  vortix  cantar; 
tha  Titeeut  dlaaipatioa  that  acta  at  a  haat  aourea  it  nora  than  offaet  by 
tha  adiabatic  axpantion  and  eoolinft  that  raaulta  from  tha  raducad  cantral 
praaaura*  Uaing  a  valua  of  v  ■  10^  cm  aac*^*  hanca 

Bl 

r  /ear  -  (a  ♦  1)  v_/a  ■  l.U  x  10^, 

ooa  finda  that  -  To  ■  0.12  dagraat  C.  Rott  ahova  that  for  the  conduct irlty 

of  air«  a  tanparatura  deficit  of  roughly  1/2  thit  magiltuda  ia  correct. 

Sinclair* a  (196U)  naaauraaanta  ahov  a  temparatura  riaa  of  from  3  to  9 
dafraaa  near  tha  center  of  duat  darlla  at  a  height  of  approximately  2  matara. 
It  thua  appa»ra  that  tha  dynatmic  reduction  of  temparatura  ia  negligible 
cceg>arad  to  the  affacte  of  tha  heated  aurfaca  boundary  layer  of  air,  vhich 
ia  praaumably  advaotad  into  tha  core  of  tha  Tortax. 
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III.  !K)P-STEADY  STATES  0?  THE  VISCOUS  VORTEX 


Vortex  ’’ith  Uni  form  Convergenee 


The  decay  of  the  vortex  vlth  unlloi’a  cocver^rence  appear*  to  have  firet 
been  obtained  bj''  Tere.iawa  {1922<  11923)  together  trlth  that  for  a  ttolel  vlth 
V  Independent  of  r  and  a  combination  "vortex**  with  v  Indapandent  of  r  In 
the  core  and  v  proportional  to  1/r  outatde  the  core. 

The  folloving  treatment  differs  only  ali^tly  from  that  of  Terasava, 
vho  vorked  in  tame  of  vortlclty  rather  than  tangential  velocltlei. 

By  setting  the  convergence  equal  to  sero  and  retaining  Sv/dt,  the 
tangential  equation  of  motion  (2.2)  becomes 


^  ♦  i  ll  .  JL. 


(3.1! 


after  transfonelng  to  the  con-dloeneional  variables 


X  ■  r(T)/r^(o),  y  »  v(t)/v_(o), 

m  Bt 

Integration  by  separation  of  variables  gives 


T  •  vt/r  ‘  (o) 

m 


y  ■  /  (x^^^x)  e'’^^dX 


vhere  e.  remains  to  be  determined. 


Prom  the  initial  conditions 
_  a  ♦  1 


(y). 


(1  - 


'T*o  a 

and  vlth  the  \iBe  of  the  Hankel  integral  transform 


0  +  1 


/  J^(X^/2^)  (1  - 


^-«x^  /2 


)  4x 


(3.2) 


(3.3) 


(3.U) 


flTlnf 
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y  “  (3.5) 

“  wo  * 

Iho  Indleatod  Intogrotioo  Moy  b«  earrioA  out  by  Boaas  of  eonfluo&t  Iqrporfoo- 

Mtrie  fbDctiona  to  flvo 


y 


a  ♦  1 
ox 


I 


(3.6) 


Thlo  raault  has  boon  plottod  for  varioua  rmliMo  of  t  Id  flgura  (^). 

’Iba  limiting  form  of  tha  solution  for  the  daei^  of  tha  Tortax  givan  bgr 
(3.6)  ms  T  bacooms  vary  Imrga  nay  ba  obtainad  amsily  fro*  sinilmrlty 
oonsidarmtioRS.*  If  tha  dacey  time  is  tekan  to  ba  grant  anou^  that  tha 
rolutioa  has  baooas  indapandaot  of  tha  initial  conditions •  diaansional 
analysis  laaus  to  a  solution  cf  tha  form. 

V  ■  vf(n)/r  vhara  n  ■  r^/vt  (3.7) 

This  allots  tha  partial  diffarantlal  aouation  (3.1)  to  ba  axprassad  as  an 
ordinary  diffarantlal  aquatXcm 


iif”(n)  +  f’<n)  ■  0 

(3.8) 

vhich  has  the  solution 

,  .  (1  . 

(3.9) 

or  in  dluensionlass  fora  o 

y .  iLiJ.  a .  F, 

V 

(3.9a) 

This  aqiiation  is  scan  to  ba  tha  asm 

a  as  tha  limit  of  aquation  (3.6)  for  t  -  * 

Tha  location  of  tha  radius  of  maxisMm  tanpantial  velocity  as  a  function 
of  tha  daoay  tiea  tMy  ba  obtainad  by  differantiating  (3.9a)  vith  raspact 
to  X  giving 


^This  appaars  to  have  baan  dona  first  by  Taylor  (1918) 


-1> 


(3.10) 


•o,  •lae«  #  •  0  «t  y  ,  It  foUow*  th«t 


♦  1  ■  • 


Tbit  bat  tha  aolutloQ  ■  2.$1  irttlcb  it  Idantieal  to  a  fron  aquation  (2.14 )* 


Tbttt 


y  .  ( 1  -  •  2-i-i  a  -  .-2*> 

a/2 

Raplaolag  a  '  by  (1  ♦  a)  and  aultlplylng  by  glrat 

V.  *  *  r^'  •  1 

Raoellingi  for  tha  constant  convarganea  aodal,  froc  (3.6) 


(3.11) 


y  .  SLIJ.  (1  .  /2(1  ♦  «^)) 

Xt  it  notad  that  tbit  aquation  hat  tha  taas  forsi  at  (3.9a)  vltb  a  ehanga 


of  Tariabla  ^ 


T'-IS-.T 


(3.12) 


ao  that  tba  loeut  of  tha  a»zina  for  (3.6)  it  alto  glrec  by  (3.11). 


"Hill  Punctioo**  Ooovar 


Iba  dacay  of  tha  vortex  {ranaratad  by  the  "hill  function"  dlatribution  of 
eoovargance  it  obtained  by  numrlcal  nathodt  free  (3.1)  utiap  at  initial  eondl- 
tiona  tba  nuaarleal  solution  of 


i 


«1»r« 


2 


•  ■  a/r.t 


KqT 

Tv 


Kaa 

2v«^ 


r/r  ,  r  ■ 


Rtvrltlnf  (3*1)  &  flnlt*  dlffartnc*  «qu^tlOB 


y  1  "*^11 

^  r_ 


2  *i  4  2h  /  “ 

X  X  x^  , 


4J 

O.i**) 


and  taking  far  boundary  eondltloaa 

y  ■  0  at  X  ■  0  for  all  t ,  and  y 


!♦!  •  * 


‘».95 


pemlta  th«  ouDaricai  aolutloa  to  b«  carried  out  until  t  becoeea  large  enougb 
that  the  almllarlty  Mclution  becoruia  ra!.id.  Stability  of  the  finite  difference 
solution  reqolres,  according  to  Hlllebranl  (I9f3)»  .P>  232»  that 


h  ^  r  ^(o) 
X  m 


(3.15) 


Xt :  V 

figure  (5)  indicates  the  results  for  an  example  vlth  ooderetely  high  eoocentra- 
tloos  of  emrergence  near  the  center •  o<  rrespondlng  to  the  point  6  ■  O.gh, 

Here  the  *^111  function**  velocity  has  been  nomalised  to  give  the  sans  value  of 
angular  scAentun  at  Infinity  as  the  non-dlncaslonal  "Burgers*  nodel."  It  Is 
noted  that  the  locus  of  the  aoximuD  tenoential  velocity  rapidly  approaches  that 
for  ths  similarity  aolutloa  as  t  becomes  large. 


P’>use*s  Model 

Bouse  (1963)  has  presented  a  vortex  model  produced  by  a  rotating 
cylinder  in  a  viscous  fluid,  by  considering  the  cylinder  to  shtink  to  an 


■15' 


infiaittly  tMLll  radius,  ifhila  slaultaasousljr  Inersaslad  la  rotational  spaad, 
to  nalntaln  a  constact  parlphsral  eiroulatloo  r^.  Aftar  aeosleretlng  tbs  rar» 
roiandlae  f^*uld  by  th«  traxMBlasloa  of  shtar  for  a  flnltt  ctneratlon  tins,  t  , 
ttM  ffansratln^  eyllodsr  Is  abrupljr  brouebt  to  rest,  and  the  surrouadinR  fluid 
allowed  to  deeelerata.  After  a  deci^  tine  of  t^  the  dlstributloa  of  teageotlal 
veloelt;*  Is  flrsa  by 

^  -  *3??  C^)  (3.1^) 

ebere 

X  •  r^Avt^,  Y  •  2*rr/rp,  T  •  (3.17) 

Here  a  differ  tat  product  ^7^  aay  be  shown  to  approach  the  Unit  1/e  as 
T  so  that  here  the  value  of  the  noxloua  circulation  becooes,  lo  the  limit, 
inversely  proportional  to  the  square  of  the  radius.* 

filnce  the  solution  of  (3<17)  satisfies  the  partial  differential  equation 
(3.1),  the  slnilarlty  solution  (3.9)  should  be  approached  as  a  linit.  Bevrltiag 
It  in  terns  of  X,  Y  and  T,  the  sisiilarity  solution  beconss 

Y  •  1  -  sxp  (^)  if  t  s  t^  (3.18) 

f 

In  this  esse  the  similarity  solution  is  approrcbed  onl:'  for  T  «  1,  that  is  for 

a  ceneratine  tins  t  very  large  in  relation  to  the  decay  tine  t.. 

8  '* 

Vertex  ^hrowth 

Roti  (1938}  hes  dlecuesed  the  growth  of  a  vortex  to  show  how  the  steady 
state  eolutlon  nes'  be  epproeehed  fron  eny  initial  dlrtributlon. 

*tiote  that  In  burgers*  oodsl,  equation  (2.13),  the  aaxinum  value  of  circulation 
renelns  et  infinite  redius. 
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St's"??  ''|7?'?»r| 


For  th»  rlscouM  rortex 

cF 

Tt 

astum*  a  folutlon  of  the  forr 

r  *  r  (o)  Mh^f  a  ■  rP(t) 
then  (3.19)  bceomee 


(3.39) 


(3.?0) 


-  I  rrr' 


(3.?1) 


vbere 


•  fr 


f  ^  and 


*‘'•37 


(3.22) 


Rearranging  (3.21)  to  give 

(P  -  I  P)  or-  -  vp3|7--  i 
and  chocs leg 

f-|r«-|p^  (3.23) 

ehere  C  la  »cm  conatant,  (3*22/  aseuBies  the  forB  of  the  steady  state  eguatloo 
vhlch  has  been  ahovn  to  bare  the  solution 

-OCr^/bv, 


r  ■  (1  -  e 

Integrating  (3.23)  gives 

-Kt 


') 


F  -  (ir  ♦ 


-1/2 


(3.2U) 


(3.25) 


where  A  la  the  constant  -of  Integral Icc. 

Pubatltutlng  r^T^  for  <^,{3*2k)  beccnee 

T  ^ 

1  -  exp  - - 


(3.26) 


as  t  '*'  •*  the  steady  state  solution.* 


•For  K  ■  0,  P-  •  (yt,(3«24y  hecones  the  ainllarlty  solutlaQ  (3.9) » 


(3.2T) 


Tor  t  ■  0  the  Inltiel  elrevdatlon  it 

fi  *  f.  |i  -  •*»  ■  r?^)] 

where  B  ■  A/C. 

Since  any  Initial  distribution  of  r  say  be  synthesised  by  the  summation  of 
solutions  (3.S7)  with  warylns  ralues  of  Bg  «3dl  initial  distributions  of  r  mist 
tend  to  the  steady  state  fontg  proirlded  k  la  constant. 

If  V  Is  not  eobstant  but  is  a  function  of  tiise,  then  (3.23)  beooiaes 


rM  p  4.  £k3  •  A 
dt  “  2  ^*2  -0 

Letticg  •  0,  then 

-  e(t)0  +00^-0 
<xx 

Nov  letting  G  «  PQ 

^loosing  P  so  that  ^  ^  ^(t)  «  0 

P  » 


and 


kSl 

Q  dt 


•  0 


dt  ■  0 


so 


-IS- 


•ad 


■  G  » 


(3.2B) 


It  !•  of  interest  to  look  at  irCt)  vhen  it  haa  the  fom  of  Ocoaoit  indf  eating 
alternate  eonvergsnce  and  divergence.  *lav  . 

-  _Binot 


where  E  is  a  constant. 

It  is  readily  seen  that  onclUatee  between  e“^  and  e,  so  that 

increases  monotonically  with  t,  and  7^  oscillates  between  decreasing 
litats  but  approaches  zero  as  t  >  so  that  it  follows  fron  (3.2U)  that 
r  -►  0  edso. 

rr.  STABILITY  0?  VOBTEX  PLOW 


the  classical  approach  of  linearized  eouatioos  for  snail  perturbations 
i 

I  in  the  mean  flow,  os  discussed  by  J.  L.  .Synge  (193B)»  C.  C.  Lin  (1933)  and 

Ictherst  nay  be  used  to  deterrine  criteria  for  hydrodynamic  stability  of  the 
^  vortex  notion.  In  the  general  case,  the  superlsQioBed  perturbation  is  aasumed  to 

9 

I  be  free  to  travel  in  any  direction  whatever  and  to  have  its  amplitude  vector 
oriented  at  random,  requiring  a  complicated  mathematical  formulation.  More 

» 

eirpl;'^,  the  perturbed  notion  may  he  considered  to  conBiat  of  eaveral  component 
traveling  waves<»nine  in  all— >of  which  three  are  coapression  vavea  end  the  oth*rs 
transverse  waves.  The  mathematical  analysis,  tentative  as  yet  owing  to  lack  cf 
cbeervatlonol  confirnation,  is  outlined  below  for  one  of  the  perturbation  moico^ 
that  of  a  sinusiodal  wave  traveling  circularly  and  having  its  amplitude  vector 
in  the  r,e  plane. 

The  motion  is  eseumed  to  consist  of  a  ase^n  component  plus  a  perturbation 

component,  so  that 

u,  ■  u  t  u,  v„  ■  V  ♦  V, 


P#  •  P  ♦  P 


(»».l) 
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vt*r«  !•  til*  total  radial  Talocity,  u  the  Man  raloe,  and  u  tba  parturbar 


A.  a _ _  _ _ ^  /a  _  -  _  a.  a —  ^  ^  M  V  ^  «b - Vi  ■  f  n'w  aa^a**  •  «V 

«ica  conpoiMati  \a  ncT^avion  inor*  OotuiJisipeDv  wj.i,a  wumt.  Ram  qwu?  ~w«u.w 


be  u  *  \i  ♦  u*,  whera  u*  if  the  perturbation  quant lt3(,but  la  not  adopted  here 
In  oi'dar  to  avoid  harlae  primed  quantltlaa  In  the  folloving  equations*)  Refer 
to  equatione  (2.1)  and  (2.2)*  where  u  is  now  rep'aced  with  u,  ■  u  ♦  u,  etc. 

The  naan  values*  u*  v  and  p*  are  seen  to  satisfy  the  ec^atlons  separately; 
and  so  these  teme  mey  be  removed  from  the  equations*  leavlns  equationa  of 
loot  ion  for  tho  perturbation.  As  it  is  assumed  that  u«u,  etc.,  the  non-ilnesr 
terms  may  be  neglected.  The  resulting  equations  for  the  rotfoa  of  small 
hcrlaontal  pe.*turbationa  ore  then 


V  8u  2*rv 

1  at 

"^2 

5  .  a  u  . 

1  au  .  1  a^u  u 

r  tF  “  r  “  p*  « 

•  ♦  w  — r*  + 

u?v  ^  V  8v  ^ 

1^.  m 

1  So 

[a^v  1  3v  1  a^v 

Zr  r  30 

eM»*  V 

pr  36  j 

[ar^  r^ae"' 

(H.2) 


(k.3) 


The  perturbed  motion  must  satisfy  the  continuity  eocdltlor 


1  fj 

a* 

r 


'!?r 


(ru)  ♦  III 


iU,k) 


Eo.uation  (U.U)  permits  a  stream  function  <'(r*6*t)  to  be  defined  such  that 
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If  equation  (t.2)  is  operated  on  by  (“ '^)  and  equatlor.  (1»,3)  b:^ 


and  the  equaticos  then  added,  the  pressure  is  eliminated  from  the  ree'utlnf 
equation,  which,  by  use  of  (4.5),  nay  be  written 


20- 


♦  iff  ^  •  V* _ ♦  ~i  ^  -  r  ♦. 


^rrt  2  ’©n  r 

T 


*rt  "  '^^mr  *  ^2  ^rr«fl  *  ^  *«ee#  *  r  Sir  *  r  '^xt« 


2v.  ▼.  V.  V 

-  ~*r«e  -  T  ♦eeo  '  T  *rr  *  7  S*  *  7  *»•’  7  *r  * 

r  r  r  r  r  r 

vtimre  th«  suffixes  Indicate  partial  deriratiTas. 

Assuming  a  solution  for  a  traTellD(^«  circular  wave. 


Iee.  ♦  lx  L. 

r*  ■ 

J 


(U.6) 


*  -  fCr).*"'*  -  (■••T) 

where  the  anpUttide  f(r)  and  the  angular  velocity  »(r)  may  be  cocplex,  but  the 
wave  DUs^er  n  is  a  real  inte/^r,  equation  (^*6)  transforms  into 


*  g-^^)  /'  ,  g  \3?ll  f-  ♦  hi  t 

r* 


(k.8) 


Which  is  strictly  appropriate  only  for  ti^  t  ■  0. 

For  purposes  of  calculation,  it  is  desirable  to  introduce  the  following 

dimensionless  quantities  in  terms  of  v  ,  the  maxiinua  mean  tangential  velocity, 

and  r  ,  the  corresponding  radius: 
m 

X  •  ^  (dlBenslooless  radius) 

*’m 

y  i  ~  (dimensionless  mean  tangectirl  velocity) 


d(x)  l.-i^f(r) 


r  V 
a  a 


(dloonsionleas  perturbation  amplitude } 


c  s 


rv_ 


I 

X 


(dineosioaless  angular  velocity  of  assn  flow) 


-2J.- 


P  E 


r  « 
n 


(din»DSiOQ7.«tt  angular  rvloelty  of  parturtetiooa) 


V  r 

B  1  -2LJE 

V 


(charactariatic  Feynolda  TIumbar  for  rortax) 

Uaing  thasa  definitional  aquation  (i).8}  nay  ba  expraaaad  in  diaenaionlaaa  fora 


X  X^  X 


InR 


Ua  -  a){r'  ♦  ♦)  - 

L  X  X 


)d 


(U.s) 


vbere  the  prlMa  nov  indicate  dlffarantiatlon  vith  raapact  to  x> 

Retumlnp  to  a  consideration  of  the  streamfunction  for  horisontal  Tortax 
flotf)  alnca  for  xean  motions  u  «  v,  a  total  streamfunction  for  t  ^sy  ba 
approximately  axpraaaed  as 


-  ?(r)  ♦  «?(r)a^^®  ■* 

.  *(r)  ♦  |i 

■  ^(r)  ■*"  V(5  ra^®  ”  (^<10) 

Cosqpariaon  of  (U.T)  and  (U.IO)  indicates  that  f(r)  ■  T(r)dri  and  slnea  £r  siay  ba 
regarded  at  an  arbitrarily  small  dispXacament ,  it  follovs  that  f(r)  «  eCr)  or* 

In  dimansionlasa  notation* 

♦  (x)  «  y(x)  ■  xa(x)  (*».ll) 

Replacing  d(x)  by  y  in  CAjustion  (U.9)  and  axprasslng  the  partuxbation  angular 
velocity  as  the  sum  of  real  and  imaginary  components •  o  •  *  ic^,  leads  to 

the  following  expressions 


•22< 


i  -  X  ♦  *  »r:  -i_  (I..U, 

o  T  y"  *  ^  2  2  2 

y'"’"'  4  ~  \  }f  .ti^,.  L  y''  4.  J.  y'  4.  y 


Bo. 


niy'*  V  y"  •  %  ^ 


(•♦,13) 

where  o^  Is  the  aniniler  pbeee  speed  of  the  perturtstions  end  Is  the  estpllfi- 
eatlon  factor  *  posit  ire  for  inltlcUy  unstchle  perturbations* 

dursers*  solution  for  the  rortex  with  unifora  horiscatsi  conrergence, 
et^'iation  (2.13) •  can  be  written  in  dimensionless  notation  as 


y -2-^(1  .  (l,.U») 

2 

whore  a  1  Kr  /cv  •  2.$1.  Curves  for  Bj.  for  various  wave  nusdsers  n,  using  this 
m  X 

expression  for  y,  are  plotted  in  figure  (6).  These  indicate  great  stability 
for  n  >  1  in  the  central  core  region,  where  tbs  nesn  flow  appiroaehes  solid 
rotation*  The  analysis  also  indicates  that  perturbations  of  wave  nuabers  1  and 
2  are  initially  unstable  for  sufficiently  Icrge  rndlli  qualitatively,  wave 
number  1  perturbations  do  appear  to  grow  in  the  outer  parts  of  vortices,  but 
quantitative  data  is  at  present  lacking. 

To  susasarize.  the  above  snalyzie  pertains  to  a  rather  special  case  of 
perturbed  motion,  where  the  only  perturbed  motion  pemitted  is  one  with 
redial  aaplltude  and  angular  oavaraant*  Cf  course,  a  truly  "laeneral**  perturba¬ 
tion  which  is  free  to  move  in  sny  direction  and  whoee  a:9litude  function 
is  throa-dlmenaional  may  have  an  entirely  different  stability  criterion, 
probcbly  highly  dependent  upon  direction  of  proposition*  Thus,  the  u\'erell 
question  of  stability  is  only  partially  aoavered  by  the  traatmrnt  above; 
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tatirwt  d«t«ndiuttion  of  vtablllty  eritcrl*  for  th»  nHmlning  aodoa  of 
porturbAtlon  >’111  romit  la  a  3  x  3  ntrir.  of  •tablllty  criteria  that  should  . 
also  describe  tbs  gecexrl  i^erturbatlon. 

V.  PLAliS  FOR  FirfflSR  STUDIES 

Three  nsla  areaues  of  study  are  iadlcatedt 

1.  Sxt«atiaa  of  the  theory  to  three-linenslonal  flow,  Incorioratlnc 
the  effects  of  conrectioo  result  lac  from  ▼arlo’is  distributions 
of  surface  beating  sad  Including  the  effecta  of  friction  In  the 

surface  boundary  layer. 

2.  Coatlnuatlcm  of  the  analysis  of  lastsbllity  for  other  perturbation 
■olea«  atteaptlng  to  ccnstruct  a  ccaplete  picture  of  dynaaic  Insta- 
b‘  ^y  asaoclated  vlth  Yarlons  vortex  models. 

3.  Developoent  of  apparatus  for  the  purpose  of  conducting  studies  of 
BO<*  s  representative  of  the  Bnallsr  scale  intense  vortices,  in 
partlciilar  dust  devils  and  fire^lnduced  whirlwinds. 
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FIG.  5!  COMPARISOH  OF  "HILL  FIBICTION"  AND  UNIFOK<  HORIZONTAL  CONVERGEHCE 
VORTEX  DECAYS. 


